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VERY LOW-MASS STELLAR AND SUBSTELLAR COMPANIONS TO SOLAR-LIKE STARS FROM MARVELS 
III: A SHORT-PERIOD BROWN DWARF CANDIDATE AROUND AN ACTIVE GOIV SUBGIANT 
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ABSTRACT 

We present an eccentric, short-period brown dwarf candidate orbiting the active, slightly evolved 
subgiant star TYC 2087-00255-1, which has effective temperature T c s — 5903 ± 42 K, surface gravity 
log(g) = 4.07 ± 0.16 (cgs), and metallicity [Fc/H] = —0.23 ± 0.07. This candidate was discovered 
using data from the first two years of the Multi-object APO Radial Velocity Exoplanets Large-area 
Survey (MARVELS), which is part of the third phase of Sloan Digital Sky Survey. From our 38 
radial velocity measurements spread over a two-year time baseline, we derive a Keplerian orbital fit 
with semi-amplitude K = 3.571 ± 0.041km s — 1 , period P = 9.0090 ± 0.0004days, and eccentricity 
e = 0.226 ± 0.011. Adopting a mass of 1.16 ± 0.11 M Q for the subgiant host star, we infer that the 
companion has a minimum mass of 40.0 ± 2.5Mj up . Assuming an edge-on orbit, the semimajor axis 
is 0.090 ± 0.003 AU. The host star is photometrically variable at the ~ 1% level with a period of 
~ 13.16 ± 0.01 days, indicating that the host star spin and companion orbit are not synchronized. 
Through adaptive optics imaging we also found a point source 643 ± 10 mas away from TYC 2087- 
00255-1, which would have a mass of O.13M if it is physically associated with TYC 2087-00255-1 and 
has the same age. Future proper motion observation should be able to resolve if this tertiary object 
is physically associated with TYC 2087-00255-1 and make TYC 2087-00255-1 a triple body system. 
Core Ca II H and K line emission indicate that the host is chromospherically active, at a level that 
is consistent with the inferred spin period and measured iVot sin j, but unusual for a subgiant of this 
T e ff . This activity could be explained by ongoing tidal spin-up of the host star by the companion. 
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1. INTRODUCTION 

Brown dwarfs (BDs) range in mass from ^13—80 
Jupiter masses and bu rn deu ter ium but not hy- 
drogen dBurrows et al.| 119971: iChabrier fc Baraffd 
[200a IBurrows et all 120011: ISpiegel et all 1201 ID. The 
first unambiguous d i scove r y of BDs (Rebolo et all 
19951: iNakaiima et all 119951: iQppenheimer et all 1199a 
Basri et all 119961: iRebolo et all 119961 ) occurred at 
the same time as the discovery of the first extra- 
solar gi ant planet orbiting a main sequence star (51 
Peg b; iMavor fc Quelol [1995T) . More than 800 BDs 
have been discovered to date (see DwarfArchives.org, 
|http: / /www.dwarfarchives.org[ ) . Most of them are 
free-floating objects and only several dozen BDs are 
companions to othe r stars (jReid fc Metchevl 120081 : 
iSahlmann et al.ll201lD . The BD desert (a paucity of BD 
companions relative to planetary or stellar companions 
within 3 AU around main-sequence FGKM stars) was 
found during high-p recision radial velocity (RV) surveys 
seeking exoplanets (jMarcv fc Butlerl 120001 ) . Since RV 
surveys are more sensitive to BDs than to exoplanets, 
this paucity is a real minimum in the mass distribution 
of close companions to solar-type stars. The California 
& Carnegie Planet Search finds a BD occurrence rate 
of 0.7 % ± 0.2% f rom their sample of ~ 1000 target 
stars (|Vogt et all 120021 : iPatel et all I2007D . and the 
McDonald Observatory Planet Search agrees, with a 
rate of 0.8% ± 0.6% fro m a search sample of 250 stars 
(jWittenmver et all 120091). To asse s s the reality of the 
BD desert. iGrether fc lineweaverl (|2006f ) performed a 
detailed investigation of the companions around nearby 
Sun-like stars. They find that approximately 16% of 
nearby Sun-like stars have close (P < 5 yr) companions 
more massive than Jupiter: 11% ± 3% are stellar, < 1% 
are BDs, and 5% ± 2% are giant planets. However, 
iGizis et all (|2001| ) suggest that BDs might not be as rare 
at w id e separations (see als o IMetchev & Hillenbrand 
12004 ). lLafreniere et all (|2007l ) obtain a 95% confidence 
interval of 1.9i®;|% for the frequency of 13— 75Mj up 
companions between 25— 250AU amongst 85 nearby 
young stars observed during the Gemini Deep Planet 
Survey. Based on an adaptive optics survey for sub- 
stellar companions. IMetchev fc Hillenbran d (2009) infer 
that the frequency of BDs in 28—1590 AU orbits around 
young solar analogs is 3.2l2'7%- 

Ostensibly, BDs are believed to form similarly to 
stars, through gravitational collapse and/or fragmen- 
tation of molecular clouds ([Padoan fc Nordl und 2004; 
iHennebelle fc Chabrierl 120081). However companions 
with mas ses up to 10 Mj„p (lAlibert "et~aTl 2005) or even 
25 M Jup (|Mordasini et al.l r2008) may form in protoplan- 
etary disks. As such, the BD desert is commonly inter- 
preted as the gap between the largest mass objects that 
can be formed in disks and the smallest mass clump that 
can collapse and/or fragment in the vicinity of a proto- 
star. In comparison, the mass function of isolated sub- 
stellar objects both in the field and in clusters appears 
to be roughl y flat in log(M) for masses down to at least 
~ 20 Mj„r, dLuhman et all 120001: IChabrierl [2002T ). Re- 
cently, lAndre et all (|2012f ) found a self-gravitating con- 
densation of gas and dust with a mass of 0.015 to 0.03M Q 
using millimeter interfcromctric observations, which sup- 
ports the idea that BDs could form the same way as stars. 



Given the occurrence rate of ~ 1% for BD compan- 
ions, a large, relatively uniform, systematic RV survey 
of a much larger sample of stars is needed to make 
further progress in understanding properties of the BD 
desert. The Multi-object APO Radi al Velocity Exop lan- 
ets large-area Survey (MARVE1S; iGe et al.l[2008h is a 
four-year RV survey of ~ 3, 300 stars with 7.6 < V < 12 
over time baselines of ~ 1.5 years per target, with a 
stated goal of < 30 m s _1 precision for the faintest 
stars. MARVE1S uses the innovative instrumental tech- 
nique of a dispersed fixed-de l ay interfe r ometer (DFDI; 
see, e.g.. lErskine fc Gel l2000l: Icl [2001 IGe et all 1200ft 
Ivan Evken et al.ll2010l : PWang et al.ll2011l . I2012allbl ) in or- 
der to simultaneously observe 60 objects at a time. By 
virtue of the large number of target stars, as well as 
the addit i on of uniform selection criteria described in 
llee et all (|2011f ). MARVE1S is well-suited to detect 
significant numb ers of rare companions. For example, 
llee et all (|201lD have recently announced MARVE1S- 
lb, a 5.9 day BD candidate around a F-type star TYC 
1240- 00945-1 located in the BD desert. 

This paper is part of a series that describes the very 
low-mass stellar and substellar compani ons to solar-like 
stars detected in the MARVE1S survey (llee et alJl2011t 
iWisniewski et al.ll2012t iFleming et al.ll2012l ). In this pa- 
per, we report a new MARVE1S BD candidate, which 
we designate MARVElS-4b, detected in orbit around the 
star TYC 2087-00255-1 (Tycho-2 star catalog: iHageTall 
[2003). In §2 we describe the observations used in this 
paper. We present stellar parameters for the star in §3 
and orbital parameters for the BD candidate in §4. We 
discuss these results and give our conclusions in §5. 

2. OBSERVATIONS AND PROCESSING 

2.1. MARVELS Radial Velocities 

TYC 2087-00255-1 was a target in the fi rst two-year 
cycle of the SDSS-III (|Eisenstein et alll20ll MARVE1S 
planet search program. This star was selected for RV 
monitoring using the presel ection methodolo gy and in- 
strumentation described in llee et ail (| 2 111 ) . The RV 
observations were taken using the Sloan Digital Sky Sur- 
vey (SPSS) 2. 5-m telescope at Apache Point Observatory 
( Gunn et aDl2006l) coupled to the MARVE1S instrument 
( Ge et all I2009T ). a 60 object, fiber-fed, dispersed fixed- 
delay interferometer. The interferometer produces two 
fringing spectra ("beams") per object, in the wavelength 
range ~ 500 — 570 nm, with resolving power R ~ 12000. 
TYC 2087-00255-1 w as observed at 23 epochs from 2009 
May 4 to 2010 July 5. llee et all (|2011f ) describe the basic 
data reduction and analysis leading to the production of 
differential RVs. The RV errors are scaled by a "quality 
factor" Q = 6.22 based on the rms errors of the other 
stars obs erved on the same SD SS-III plate as TYC 2087- 
00255-1 (jFleming et al.ll2010D . The differential RV mea- 
surements for TYC 2087-00255-1 from MARVE1S are 
summarized in Table [TJ Note that a constant velocity 
term has been subtracted to account for the instrument 
offset (see ij4.1l for more detail). 

2.2. Spectra For Stellar Characterization 

In pursuit of precise stellar parameters for the pri- 
mary, optical (~ 3500 - 9000 A) spectra of TYC 2087- 
00255-1 were obtained using the FEROS high resolution 
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TABLE 1 

SUMMARY of MARVELS RADIAL VELOCITIES 
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{R = 48000) spectrograph (jKaufer et al.lll999| ) mounted 
on the MPG/ESO 2.2-m telescope in La Silla on 2010 Au- 
gust 2. FEROS spectra were analyzed using the online 
FEROS Data Reduction System (DRS) and the standard 
calibration plan, where bias, flat-field and wavelength 
calibration lamp frames are observed in the afternoon. 
Three 2400s exposures were combined to yield a S/N of 
~ 230 per one-dimensional extracted pixel at 6600 A. 

Additional spectroscopic observations around the Ha 
line were also secured with the coude spectrograph of 
the 1.60m telescope at Observatorio do Pico dos Dias, 
Brazil on 2010 August 17. The resolution was set to 
R =18 000 and the S/N per pixel was ~ 80. Data re- 
duction was carried out by the standard procedure using 
IRAF. After usual bias and flat-field correction, we sub- 
tracted the background and scattered light and extracted 
one-dimensional spectra. No fringing was present in our 
spectra. The OPD coude single order spectrograph in- 
troduces no necessity for blaze function corrections and 
thus the line profile is easily normalized, lending itself to 
accurate analysis of the temperature profile of the stellar 
atmosphere by fitting the observed profile to theoretical 
calculations. This will be used below to infer an inde- 
pendent T e g estimate as well as the evaluation of the 
chromospheric radiative losses in the Ha line core. 

2.3. Additional Radial Velocity Observations 

High-resolution spectra were collected with the Spec- 
trograf o di Alta Reso l uzion e Galileo (SARG) spectro- 
graph (jGratton et al.l 120011 ) at the 3.58m Telescopio 
Nazionale Galileo (TNG) for additional RV measure- 
ments from 2010 August 27 to 2011 August 17. This 
spectrograph provides R ~ 57, 000 spectra spanning a 
wavelength range of 462—792 nm. The spectra were 
reduced using standard IRAF Echelle reduction pack- 
ages. Frames were trimmed, bias subtracted, flat-field 



corrected, aperture-traced and extracted. We obtained 
15 epochs of observations with an iodine cell and an addi- 
tional epoch without the iodine cell to serve as a stellar 
template. The exposure time for each epoch is 20 to 
30 minutes. The signal-to- noise ratio (S/N) per resolu- 
tion element at 550 nm is ~ 60 — 250, and one resolu- 
tion element is sampled by 4.9 pixels. A total of 15 RV 
data points were derive d using the iodine cell technique 
(|Marcv fc Butler! f2000h . Each of 21 SARG orders be- 
tween 504 and 611 nm was subdivided into 10 sections, 
and then RVs were measured from these components. 
Following a 2 — a clip, the measurements were averaged 
to produce the RVs. The resultant differential RVs are 
summarized in Table 2. Note that a constant velocity 
term has been subtracted to account for the instrument 
offset (see gj} 

TABLE 2 

SUMMARY of SARG RADIAL VELOCITIES 



HJD 


RV 
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0.032 



2.4. Diffraction-limited Imaging 

High angular resolution lucky images (Lis, observa- 
tions taken at very high cadence to achieve nearly- 
diffraction-limited images fro m a subsample of th e total) 
were obtained using FastCam ()Oscoz et al. II2008T) on the 
1.5 m TCS telescope at Observatorio del Teide, Spain. 
The primary goal of these observations was to search for 
companions at large separations that could pollute spec- 
troscopic RV observations of the targets. The LI frames 
were acquired on 2011 May 8 and 2011 July 1 in the I 
band, covering ~ 21" x 21" on the sky. A total of 60000 
images, each corresponding to 60 ms integrations, were 
taken on 2011 May 8, and a total of 64000 frames, each 
corresponding to 50 ms integrations, were taken on 2011 
July 1. 

To further assess the multiplicity of TYC 2087-00255- 
1, we acquired AO images using NIRC2 (instrument 
built by Keith Matthews) on the Keck II telescope 
on 2012 August 25 UT. TYC 2087-00255-1 is bright 
(V=10.6) and served as its own on-axis natural guide 
star. NIRC2 is a high-resolution near-infrared camera 
that provides a plate scale (whe n operating in na rrow 
mode) of 9.963 ± 0.006 mas/pix (|Ghez et al.l (2008h and 
10.2 x 10.2 field of view. Our observations consist of 
dithered frames taken with the K' band (A c = 2.12/im, 
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Fig. 1. — Continuum-normalized, high resolution FEROS spec- 
tra of TYC 2087-00255-1, which is a GOIV star. For comparison 
purpose, also shown in this plot are t he spectrum of a GOIV star 
HD 098630 from Elodie spectra library (Prugnicl & Soubiran 2001) 
and the synthetic spectrum of TYC 2087-00255-1 calculated us- 
ing SME package (please see the text for more information). The 
FEROS spectrum and the synthetic spectrum have been shifted in 
the y-axis direction for display purpose. 

AA = 0.35/im) filter. The total on-source integration 
time was 47.5 seconds. 

2.5. SuperWASP Photometric Data 

To check for intrinsic photometric variability indicat- 
ing stellar activity and search for possible transits of the 
companion, we extracted photometric time series data of 
TYC 2087-00255-1 f rom the SuperWASP public archive 
(|Butters et all 120101) . The WASP instruments provide 
flux measurements for millions of stars using wide-angle 
images of the night sky over a band pass of 400-700 
nm defined by a broad-band filter. Eight cameras on 
each instrument provide images covering approximately 
7.8 x 7.8 degrees using Canon 200mm f/1.8 camera lenses 
and e2v 2048 x 2048 CCDs. Synthetic aperture photom- 
etry using an aperture radius of 49 arcseconds at the 
position of targeted stars is performed on the images 
(jPollacco et alll2006h . A total of 18935 aperture photom- 
etry data points, each taken with a 30 second integration 
time, were available from the SuperWASP public archive. 
These data were taken between 2004 May 2 and 2008 Au- 
gust 10. Systematic errors caused by spatially localized 
flat-fielding, errors in the vignetting correction near the 
edge of the field of view, bright moonlight contamina- 
tion, bad weather and other as-y et unidentified reasons 
do ex ist in SuperWASP data sets ((Collier Cameron et al.l 
2006). There are large systematic errors in some of the 
data sets, and we choose to use data points with rela- 
tive flux errors smaller than 0.01 in our further analysis, 
to limit systematics. The final selected data sets have a 
total of 11932 photometric data points. 

3. TYC 2087-00255-1: THE STAR 

3.1. Stellar Parameters 

We have used two independent pipelines (referred to 
as "BPG" and "IAC") to derive the stellar parameters 
from the high resolution FEROS spectra. Both methods 
involve analysis of the equivalent widths of Fe I and Fe II 
lines to balance the excitation and ionization equilibria 
of these features. These two me thods are described in 
detail in [Wisnicws ki et all (|2012j ) . Both pipelines pro- 
duce values of effective temperature T e g, surface gravity 




3931 3932 3933 3934 3935 3936 3937 

wavelength (A) 



Fig. 2. — Continuum-normalized, high resolution FEROS spec- 
tra of TYC 2087-00255-1 centered on the Ca II K line. The excess 
emission from the Ca II K line core indicates that TYC 2087-00255- 
1 is a chromospherically active star. Also shown on this plot are 
several chromospherically quite stars with similar st ellar parame- 
ters to TYC 2087-00255-1 for comparison purpose (G hezzi et alj 
l2010al lbT). 

log(g) and metallicity [Fe/H] mutually consistent. Be- 
cause both pipeline determinations are mutually consis- 
tent, we average these two sets weighted by their own 
pipeline errors to determine the final stellar parameters, 
shown as 'combined results' in Table [3] For each stel- 
lar parameter, we add in quadrature a systematic er- 
ror of 18 K, 0.08 and 0.03 for T cS , log(g) and [Fe/H], 
respectively, in addition to the inte rnal errors inher- 
ent f rom the two pipeline results (see [Wisnicwsk i et al.1 
120121 ). The results are summarized in Table [3] The 
adopted stellar parameters are T ff = 5903 ± 42 K, 
log(g) = 4.07 ± 0.16 (cgs), and [Fe/H] = -0.23 ± 0.07. 
So TYC 2875-00255-1 is a GOIV type star. In Fig. □ we 
show the continuum-normalized high resolution FEROS 
spectra of TYC 2087-00255-1 together with the spectra 
of HD 0986 30, a known GOIV type star from Elodie spec- 
tra library (jPrugniel fe Soubiran! I200"TI ). For comparison 
purpose, we also show synthetic spectra of TYC 2087- 
00255-1 in Fig. [T] The synthetic spectra is calculated 
using software package Sp ectroscopy Made Easy (SME; 
iValenti fe Piskunovl 119961 ) and the stellar parameters 
from the above analysis. 

Continuum-normalized, high-resolution FEROS spec- 
tra of TYC 2087-00255-1 centered on the Ca II K line 
is shown in Fig. [3] Also shown in this plot are three 
stars with similar stellarparameters as TYC 2087-00255- 
1 jGh czzi etaLl EoiOiB)- The excess emission from 
the Ca II K line core indicates that TYC 2087-00255- 
1 is more active than stars with similar stellar parame- 
ters. Using the high- resolution FEROS sp ectra and the 
method described in Uenkins et al.1 (|2008f) , the chromo- 
spheric Ca II HK activity index of TYC 2087-00255-1 is 
logi? HK = —4.58, with a calibration rms error of 0.03. 

The Ha profile of TYC 2087-00255-1 is shown in Fig. [3] 
superimposed over the solar one, the latter obtained as 
a disk-integrated spectrum from Ganymede (which re- 
flects light from the Sun) under the same observational 
conditions. The shallower wing profile is apparent, trans- 
lating into a lower T cff for TYC 2087-00255-1, as is 
the much stronger line core filling, interpreted as ad- 
ditional chromospheric fill-in. Note that the Ha line 
core is substantially broader in TYC 2087-00255-1, inter- 
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TABLE 3 

Spectroscopic Parameters of the Star TYC 2087-00255-1 



T e g 


log(s) 


[Fe/H] 


it 


Notes 


(K) 


(cgs) 




(km s 1 ) 




5805 ± 71 


4.02 ± 0.18 


-0.24 ± 0.10 


1.74 ± 0.10 


BPG (ESO 2.2-m) 


5941 ± 44 


4.15 ± 0.22 


-0.23 ± 0.08 


1.687 ± 0.055 


IAC (ESO 2.2-m) 


5903 ± 42 


4.07 ± 0.16 


-0.23 ± 0.07 


1.70± 0.06 


combined results 



TABLE 4 

Derived Parameters of the Star TYC 2087-00255-1 



Parameter Value 
Spectral Type GO IV 

Mass 1.16 ± 0.11 M Q 

Radius 1.64±O.37i? 
Age 5.5 Gyr 



A V 0.12+-- 
Distance 218 ± 14 pc 
log<K -4.58 



preted as yet another con firmation of the subgiant sta- 
tus of TYC 2087-00255-1 (iPasauini fc Pallaviciml HMUt 
ILvra fc Porto de~M cllo 2005^ The chromospheric loss 
in the Ha core of TYC 2 087-00255-1 was also evaluate d 
under the prescription of ILvra fc Porto de Mellol (|2005[) . 
using as input the spectroscopic atmospheric parameters. 
We have derived a total chromospheric flux of 13.4 x 10 s 
erg cm" 2 s" 1 inside the Ha line core, the estimated er- 
ror in the ILvra fc Porto de Mello (2005|) procedure being 
~ 0.5 x 10 5 erg cm" 2 s _1 . This value is probably also 
contaminated by veiling from the companion, and should 
be taken as an upper limit. Nonetheless, it is over four 
times the expected fl ux from typical subgiants in this 
T cff range (Fig. 3 of ILvra fc Porto de Mellol 12001 and 
lies between the average chromospheric radiative losses 
of Pleiades (age ~ 100 Myr) and Ursa Major group stars 
(~ 400 Myr). This flux is in very good agreement with 
the Ca II K line profile (which being in the UV is proba- 
bly free from the secondary's contamination) in pointing 
to a very high activity level in TYC 2087-00255-1, com- 
patible with a solar-type stars no more than a few hun - 
dred million years of age (|Lvra fc Porto de Mellol [20051 ) . 
Two independent spectroscopic chromospheric indicators 
therefore confirm TYC 2087-00255-1 as much more ac- 
tive than expected from its subgiant status. 

The spectral energy distribution (SED) was con- 
structed for TYC 2087-00255-1 in Fig. H using near 
UV (GALEX.lMorrissev et al J 12001 . optical (H0 g et al.l 
20001 IKharchenko fc Roe serl 120091). near IR (2MASS, 
Cutri et all I2003T) and IR (WISE, IWright et all fioWk . 
Cutri et all 120121 ) photometric data. These photomet- 
ric data are presented in Table [5J The data were 
fit with fluxes from a NextGen model atmosphere 
(|Hauschildt et al.lll999l ). We limited the maximum line- 
of-sight extinction to be Ay < 0.15 base d on the analy- 
sis of dust maps by iSchlegel et al.l ([1998) . The resultant 
parameters, T cS = 5700 ± 200 K, log(g) = 4.0 ± 0.5, 
[Fe/H] = -0.5 ± 0.5 and A v = 0.12±g;g|, agree within 1- 
er of the results found via analysis of our high resolution 
spectroscopy. In the above analysis, we did not constrain 
any of the fit parameters except for Ay- We did another 
fit where we forced T c g, log(g) and [Fe/H] to the spec- 
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Fig. 3.— Continuum-normalized Ha profile of TYC 2087-00255- 
1, with a spectroscopic resolution of 18 000. Also shown on the 
plot is Ha profile of the Sun for comparison. The solar spectrum is 
measured from the reflected light of the Sun by Ganymede, taken 
together with TYC 2087-00255-1 using the same instrument. The 
shallower wing profile and the stronger core filling in the Ha line 
when compared to the Sun suggests an effective temperature lower 
than the spectroscopic temperature, which is also found by the 
SED analysis. 
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Fig. 4. — The observed SED from the near-UV through the IR 
for TYC 2087-00255-1, along with a best-fit NextGen model at- 
mosphere. Blue points represent the expected fluxes in each band 
based on the model, red horizontal bars are the approximate band- 
pass widths, and red vertical bars are the flux uncertainties. There 
is potentially some GALEX FUV excess indicating this star is chro- 
mospheric active. The resultant fundamental stellar parameters 
from this fit agree with the stellar parameters determined from the 
stellar spectra to within 1-cr. 

troscopically determined values, which provide a more 
robust estimate of Ay = 0.12 ± 0.03. Using this total 
extinction estimate, and adopti ng a V ban d bolometric 
correction BC V of -0.19 ± 0.02 (|Coxll2000D . we estimate 
the distance to this system is 218 ± 14 pc. 
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TABLE 5 

Photometric Measurements of the Star TYC 2087-00255-1 



Parameter 



Value 



Note 



galFUV 

galNUV 

B 

BT 

V 

VT 

J2M 

H2M 

K2M 

WISE1 

WISE2 

WISE3 

WISEA 



21.285±0.506 
15.511±0.013 
11.203±0.059 
11.32± 0.06 
10.553 ± 0.048 
10.58 ± 0.04 
9.29 ± 0.02 
8.96 ± 0.03 
8.88 ± 0.02 
8.839 ± 0.024 
8.864 ± 0.022 
8.792 ± 0.027 
9.169 ± 0.454 



GALEX (Morrissey et al. 2007) 
GALEX (Morrissey et al. 2007) 
Kharchcnko & Rocscr 
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Fig. 5. — The evolutionary track for an object with M = 
1.16M Q , at [Fe/H] = -0.23. Ages of 1, 2, 5, 5.5, and 6Gyr are 
indicated as dots. The possible tracks for ±1ct deviation in the 
mass are shown by the shaded region. The stellar parameters for 
TYC 2087-00255-1, with l-a error bars, are shown by the cross. 

3.2. Stellar Mass and Radius 

We determine the mass and radius of the parent star, 
TYC 2087-00255-1, from T off , log( q), and [Fe/H] usin g 
the empirical polynomial relations of lTorres et al.l (12010), 
which were derived from a sample of eclipsing binaries 
with precisely measured masses and radii. We estimate 
the uncertainties in M* and i?„ by propagating the uncer- 
tainties in T c ff , log(g), and [Fe/H] (see T able 111) using the 
covariance matrices of the iTorres et al.l (|2010f) relations 
(kindly pro vided by G. Torres). Since the polynomial 
relations of lTorres et "all (|2010l) were derived empirically 
the relations were subject to some intrinsic scatter, which 
we add in quadrature to the uncertainties propagated 
from the stellar par ameter measureme nts (eri ogm = 0.027 
and cr logr = 0.014; rTbrres et al.ll20Tol) . The final stellar 
mass and radius values obtained are M* = 1.16±0.11 M© 
and i?» = 1.64 ± 0.37 R Q (see Table |4}. 

3.3. Evolutionary State 

In Fig. [S] we compare the spectroscopically measured 
T cff and log(g) of TYC 2087-00255-1 (red error bars) 
against a theoretical stellar evol utionary track from th e 
Yonsei-Yale ("Y 2 ") model grid (|Demarque et al.ll2004D . 
The solid curve represents the evolution of a single star 
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Fig. 6.— Top: Phase-folded light curve for TYC 2087-00255- 
1 at a period of 13.16 day from SuperWASP. Bottom: Lomb- 
Scargle periodogram of the SuperWASP data, showing no evidence 
for any significant periodicities around P = 9 days (frequency of 
0.11 days -1 , the orbital period of MARVELS-4b). Instead there 
is evidence of a significant period at 13.16 days (frequency of 0.076 
days -1 ), which is likely the rotational period of TYC 2087-00255-1 
as tracked by rotational modulation of starspots. There are several 
peaks around the 1-day period, which are daily aliases. 

of mass 1.16 ± 0.11 A/© and metallicity of [Fe/H] = 
—0.23 ± 0.07. The dashed curves represent the same 
evolutionary track but for masses ±0.08 M©, which rep- 
resents the 1-tr uncertainty in our derived mass. The 
filled gray region between the mass tracks represents the 
expected location of a star of TYC 2087-00255-1's mass 
and metallicity as it evolves off the main sequence. The 
spectroscopically measured T e g, log(<?), and [Fe/H] place 
TYC 2087-00255-1 in the subgiant phase, prior to the 
base of the red giant branch, with an estimated age of 
- 5.5 Gyr. 

3.4. Stellar Rotation Period and Rotational Velocity 

In this section, we will use three different ways to es- 
timate the rotation period of TYC 2087-00255-1. First, 
we find a sinusoidal variation in the 2004 SuperWASP 
photometry data with a period P = 13.16 ± 0.01 days 
and amplitude of 9 mmag. In Fig. [5] we have shown a 
phase-folded plot of the 2004 SuperWASP data. Since 
the SuperWASP photometry data used an aperture ra- 
dius of 49 arcseconds, we used the SIMBAD website 
to check for bright stars inside this 49 arcseconds aper- 
ture which could account for this 9 mmag variation and 
found several stars with V mag > 16.2. Such faint stars 
could not produce such a 9 mmag variation around the 
V mag = 10.6 TYC 2087-00255-1 unless their luminosities 
varied by 100%, which is quite unlikely. We can explain 
this 13.16 d period as the rotation period of the host star 
and the 9 mmag sinusoidal variation by the rotational 
modulation of star spots. 

The second method to derive the rotation period 
is to use the chromospheric a ctivity-rotation relation 
(jMamaiek fc Hillenbrand! 12008). The chromospheric 
Ca II HK activity index of TYC 2087-00255-1 derived 
previously is logi? HK = —4.58. The corresponding rota- 
tion period is 11.8 d, with an estimated error of 2.4 d 
from this calibrated relation. This rotation period agrees 
with the one derived above (13.16 d) using the Super- 
WASP photometry data. However, we should note here 
that this chromospheric activity-rotation relation is de- 
rived for solar-type dwarf stars. Since TYC 2087-00255-1 



7 



■ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


I 1 ' 


ii iii hi ini iti 1 1 1 in 1 1 1 1 iii 

II I P || 1 




1 1 III 1 1 1 III 


1 1 III 1 1 1 III 1 1 1 [If 

| |l j 


I 'l 1 




Vr<- 


r i 
- 1 ( 




1 



5000 5200 5400 5600 5800 

HJD-2450000 



Fig. 7.— Top: Keplerian orbital solution for TYC 2087-00255- 
1. Open triangles are MARVELS discovery data, open squares are 
TNG data. Bottom: The residuals between the data points and 
the orbital solution. 
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Fig. 8. — Top: Phase-folded Keplerian orbital solution and RV 
residuals for TYC 2087-00255-1. Open triangles are MARVELS 
discovery data, open squares are TNG data. Bottom: The resid- 
uals between the data points and the orbital solution. 

is a slightly evolved subgiant, it may not be appropriate 
to use this relation here. 

Thirdly, we use the equation 2nR*/v ro t s'mi to esti- 
mate the rotation period of the host star. We measured 
the projected rotational velocity iwsini of TYC 2087- 
00255-1 by comparing our high resolution SARG and 
FEROS spectrum to broadened versions of Kurucz AT- 
LAS synthetic spectra. We used the atmospheric pa- 
rameters derived in §3.1 and fixed the macro turbulence 
velocity to V m acro = 4.2 km s" 1 based on Equation 
(1) from I Valenti fc Fischerl ([20051 ) . We find v rot smi = 
9.2 ±2.0 km s _1 using the observed SARG spectrum and 
v ro t sini = 10.1 ±0.9 km s -1 using the observed FEROS 
spectrum. Using the derived radius of TYC 2087-00255- 
1 -R* = 1.64±0.37i? Q , the corresponding rotation period 
are 9.1±2.9 d and 8.3±2.0 d. These are about 1.4-sigma 
and 2.4-sigma from the 13.16 day period inferred from 
photometric data, suggesting that the star is likely close 
to edge on, as a smaller inclination would imply shorter 
periods and thus larger discrepancy with the photometric 
period. The slight tension even assuming sini = 1 may 
arise from systematic errors in the estimate of tVotsinz, 
including an incorrect assumed value for V ma cro- 

4. TYC 2087-00255-1: THE COMPANION 
4.1. Keplerian Orbital Solution 



Radial velocities derived from MARVELS and SARG 
data were used to fit Keplerian orbital parameters. Since 
there are star spot features seen from the light curve of 
TYC 2087-00255-1 (see we added an additional 

stellar 'jitter' ter m (gut t er) in our Keplerian orbital model 
as suggested by iFordl ()2006[ ) to account for any addi- 
tional noise induced by the stellar activity. We performed 
two Keplerian orbital fits using th e Markov C hain Monte 
Carlo method (MCMC, see, e.g.. iFbrdl 120061 ) . First, we 
use only the MARVELS data to fit the Keplerian orbit, 
then we combine the MARVELS and SARG RV data to 
do a joint fit. For the details of our o ne planet/BD RV 
model, please see §2 of iGregorvl ([2007D . The best-fit pa- 
rameters from the two fits are presented in Table [6] and 
they agree with each other quite well. 

As expected, the fit of the combined MARVELS and 
SARG RV data has relatively smaller uncertainties com- 
pared to the fit with just the MARVELS data, there- 
fore we use this combined fit as our final Keplerian or- 
bital solution. In this fit, a constant systematic veloc- 
ity term is included for each of the two instruments 
to account for the offset between the observed differ- 
ential RV data and the zero-point of the Keplerian 
RV model (-0.803 ± 0.035 km s" 1 for MARVELS and 
-353 ± 0.029 km s" 1 for SARG). The RV data shown 
in Table Q] and Table [2] are the RVs after subtrac- 
tion of these two constant systematic velocity terms. 
The final Keplerian orbit of MARVELS-4b has a pe- 
riod P = 9.0090 ± 0.0004 days, e = 0.226 ± 0.011 and 
semi-amplitude K = 3.571 ± 0.041km s" 1 . This solution 
is shown in Figs [7] and |8] together with the MARVELS 
and SARG RV data. The residuals shown in these two 
plots could not be explained only by the errors in our 
RV data. A stellar jitter term cotter = H2 ± 17m s _1 
is required in our fitting to explain these residuals. In 
our previous p aper ([Fleming et al.|[20l0l: iLee et alJ[20lH 
iWisniewski et al.ll2012t ).' we did not include a stellar 'jit- 
ter' term in our Keplerian fitting. Instead we re-scaled 
the error bar of MARVELS RV data to force the reduced 
chi-square to be 1. Since MARVELS-4b is orbiting an 
active star which has spots activity, we choose to include 
this 'jitter' term and did not try to rescale the MAR- 
VELS and SARG RV error bar. We note here that both 
methods (including a 'jitter' term or rescaling the RV 
error bar for the reduced chi-square to be 1) yield Ke- 
plerian orbital parameters consistent with each other in 
the la range. We will discuss more about this stellar 
jitter in 35.31 Using the derived value of M* in 33. 2[ 
we estimate a minimum mass (i.e., for sin a = 1 where a 
is the line-of-sight orbital inclination) for the companion, 
MARVELS-4b, as m min = 40.0±2.5 M Jup , where the un- 
certainty is dominated by the uncertainty in the primary 
mass. Assuming the system is edge on, the semimajor 
axis of this system is a = 0.090 ± 0.003 AU. 

4.2. Search for Possible Transit Signals 

Based on the orbital parameters of the companion, 
we have calculated the a priori tr ansit probability to 
be 7.8 % using Equation (5) from iKane fc von Braunl 
(2008). The expected duration of a central transit is 
i~ R*P/(ira) = 2.94 hr, and the expected depth is 
- (R/R*) 2 = 0.004(i?/ J R Jup ) 2 , where R is the radius of 
the companion. We phase folded the SuperWASP data 



TABLE 6 

MARVELS-4B: Parameters of the Companion 



Parameter 


MARVELS+TNG 


MARVELS 


Minimum Mass 


40.0± 2.5 M Jup 


39.9± 2.5 M Jup 


a 


0.090 ± 0.003 AU 


0.090 ± 0.003 AU 


K (km s" 1 ) 


3.571 ± 0.041 km s" 1 


3.563± 0.073 km s" 1 


P(d) 


9.0090 ± 0.0004 d 


9.0105 ± 0.0024 d 


e 


0.226± 0.011 


0.233± 0.022 


w 


4.086± 0.041 


4.077±0.081 


""jitter 


Q.112± 0.017 km s" 1 


0.152± 0.046 km s" 1 


^prediction for transit (HJDtjxc) 


2455549.629 ± 0.056 


2455549.715 ± 0.082 


Tperiastron (HJDtjxc) 


2455552. 797±0. 083 


2455552.851±0.147 
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Fig. 9. — Composite images showing the results of different LI 
thresholds applied to frames acquired with FastCam on May 8th 
2011 and July 1st 2011. Each panel covers ~ 5.5 X 5.5 square 
arcseconds, centered on TYC 2087-00255-1. To ease visualization, 
each image has been normalized to a peak value of unity. This set of 
images illustrates the gain in angular resolution close to the target 
location when applying highly restrictive LI thresholds, but at the 
cost of lowering the contrast achieved at larger angular distances 
from target location. 

and searched for a transit signal at the expected transit 
time, but no evidence for a transit with a period ~ 9 days 
was found. We could rule out existence of a transit sig- 
nal with depth greater than 0.014 mag at 3-cr confidence 
level. Our conclusion is that the SuperWASP data is 
not sensitive enough to detect a 4 mmag transit signal. 
Further photometric follow-up observations are needed 
to rule out or confirm this transit signal. 

4.3. Search for Possible Stellar Companions Using 
Lucky Imaging 

We use the lucky images to search for, and place con- 
straints on, any possible undetected stellar companions 
at large separations. The data were processed using a 



custom IDL software pipeline. After identifying cor- 
rupted frames due to cosmic rays, electronic glitches, etc., 
the remaining frames are bias corrected and flat fielded. 
Lucky image selection is applied using a variety of selec- 
tion thresholds (best X%) based on the brightest pixel 
(BP) method. The selected BP must be below a spec- 
ified brightness threshold to avoid selecting cosmic rays 
or other non-speckle features. As a further check, the 
BP must be consistent with the expected energy distri- 
bution from a diffraction speckle under the assumption 
of a diffraction-limited PSF. The frames are sorted from 
brightest to faintest according to the brightness of their 
BPs and the brightest X% are then shifted and added to 
generate a final image. In the different panels of Fig. [9j 
we show the composite lucky images generated from con- 
sidering only the best {1, 5, 15, 25, 50, 80}% of the image 
frames. Each panel covers ~ 5.5 x 5.5 arcsec 2 centered on 
TYC 2087-00255-1. Restricting the images used to the 
best 1% improves the angular resolution with respect to 
the loose LI selection (i.e. the top 80% of the images 
image), but at the cost of increased noise at larger sep- 
arati ons from the target. F ollowing the same procedure 
as in iFemema et al. I (|2011[) , we compute the 3-cr "best 
LI detectability curves" for each of the two nights' data, 
which are shown in the left panel of Fig. [101 We use 
the following procedure to place a 3-cr upper limit on the 
mass of a possible undetected stellar companion: 

1. From the SED fitting of TYC 2087-00255-1 in ^3~T1 
that provides an estimate of the extinction and dis- 
tance, we estimate the absolute / band magnitude 
of TYC 2087-00255-1 to be M 7 = 2.74. 

2. Knowing Mj and the 3-cr detectability curves (see 
the left panel of Fig. [TOT) allows us to construct 
the Mi versus p (angular distance from TYC 2087- 
00255-1) curve for TYC 2087-00255-1. This curve 
provides the upper limit of absolute / band mag- 
nitude for any undetected stellar companion at the 
3-cr level. 

3. Although TYC 2087-00255-1 is identified as a sub- 
giant, the fact that it is not very massive makes it 
plausible that any stellar companion with the same 
age and smaller mass will still be a main sequence 
object. These facts j ustify the u s e of th e conversion 
from I to V band in iMamaiekl (|2010f ). 

4. From the My versus p curves we employ the empir- 
ical Mass-Lu minosity relationships (MLRs) from 
the literature (jHenrv et al. 1 119991 : iDelfosse et al. I 
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Fig. 10. — Left: Best 3cr I-band detectability curves from lucky 
images achieved on 2011 July 1 and 2011 May 8. We see the 
quality of both nights are comparable. Right: Conversion of 3-<r 
detectability curves into mass sensitivities using empirical Mass- 
Luminosity relationships in the literature. See the main text for 
details. 

20001; IHenrv I l200i IXia et al. I [20081 : IXia fe Fu I 



2010) to derive the 3-cr upper mass limit for any 



undetected companion as a function of angular dis- 
tance to TYC 2087-00255-1. 

The results of applying the above procedure to our LI 
data are shown in the right panel of Fig. [TOl 

4.4. Search for Possible Stellar Companions Using A 

Imaging 

In this section we use the acquired Keck AO images to 
search for possible stellar companions around TYC 2087- 
00255-1. The AO images were processed by replacing hot 
pixel values, flat-fielding the array, and subtracting ther- 
mal background noise. No companions were identified in 
individual raw frames during the observations. However, 
upon stacking the images we noticed a point source to the 
north-east of TYC 2087-0025501. Fig. HT]shows the final 
processed K' image. The candidate is 6.49 mags fainter 
than the primary star in K'. We measure an angular sepa- 
ration of 643±10 mas and posit ion angle 2 7.1 ±0.1 de gree. 
Assuming an age of 5 Gyr, the Bara ffe et all (,19981 ) the- 
oretical evolutionary models predict a mass of O.13M if 
the candidate is physically associated at a distance of 218 
pc. With a proper motion of (/ i a cos S, us) = ( —2.9, 39.8) 
mas/yr for TYC 2087-00255-1 ([Hog et al.ll2000h . it will be 
possible to determine whether this candidate is a tertiary 
companion in less than one year with NIRC2. 

5. DISCUSSION 

5.1. Does MARVELS~4b Reside in the Low-mass Tail of 
the Stellar Formation Process? 

The mass of BD overlaps both with low-mass stars 
formed from collapse/fragmentation of molecular cores 
and massive planetary companions formed in protoplan- 
etary disks. It is still not clear which mechanism domi- 
nates the formation of BDs. By extrapolating the com- 
panion mass functio n from both the exopla n et sid e and 
low-mass star side, iGrether fe Lineweaverl |2006) find 
the minimum number of companions per unit interval 
in log mass is 3ll^gMj up . iSahlmann et al.1 ([201 ll ) see 
evidence for a bimodal distribution in BD masses, with 
the gap between 25 and 45 Mj up almost entirely de- 
void of objects. They suggest that the less-massive 



Fig. 11.— Keck K' band AO image of TYC 2087-00255-1. A 
candidate 6.49 mags fainter than the primary star in K' band is 
identified. It has an angular separation of 643±10 mas and position 
angle 27.1 ±0.1 degree (pointed by an arrow in the plot). 

group may represent the high-mass tail of the plane- 
tary distribution. If true, the maximum mass of giant 
"planets" should be around 30Afj up . The BD candidate 
MARVELS-4b reported in this paper has a minimum 
mass of 40.0±2.5M Jup , which suggests that MARVELS- 
4b more likely formed like stars through collapse and/or 
fragmentation of molecular cores. 

5.2. Activity, Rapid Rotation and Stellar Spin Evolution 

TYC 2087-00255-1 is in its subgiant phase with an 
estimated age of ~ 5.5 Gyr and chromospheric Ca II 
HK in dex logi? HK = —4.58 (see 33. ip . Uenkins et al.l 
(|2011|) have studied chromospheric activity indices for 
more than 850 FGK-type dwarfs and subgiant stars, and 
find the distribution of activity indices (log i? HK ) for their 
subgiant sample can be fit by a Gaussian centered around 
-5.14, with cr = 0.06. The activity index for TYC 2087- 
00255-1 is log i? HK = —4.58, which makes it an unusually 
active subgiant star (see also § 13. If) . 

For cool stars with T c g < 6500 K, chromospheric ac- 
tivity is generated through a stellar magnetic dynamo, 
which is related to the rotational velocity and rotation 
period of the star. The star's rotation period increases 
with age throu gh mass loss in a magnetized wind ( "mag- 
netic brakin g " iSchatzmanl Il962t IWeber fe Davis] 119671; 
Mestell 119681 iSkumanichl Il972t lEpstein fc Pinsonneaultl 



20121 ) . and as such, its chromospheric activity level 
is expected to also de cay with age, a phenomenon 
that has been observed (| Wilson! Tl963t ISkumanichl 119721 : 
iSoderblom et alj I1991D . For an evolved star such as 
TYC 2087-00255-1, the rotation period should be rela- 
tively large, >30 days due to magnetic braking and con- 
servation of angu l ar momentum as the radius e xpands 
(|Skumanichlll972l : lEpstein fe Pinsonneaultll2012l ). How- 
ever, we find the rotation period (13.16 days) is much 
shorter than this value. One possible explanation is that 
the tidal interaction with the companion has spun up 
the star and keeps the stellar magnetic dynamo active. 
In this scenario, the tidal interaction transfers orbital 
angular momentum into stellar rotational angular mo- 
mentum. We will explore the coupled roles of radial ex- 
pansion and tidal evolution, and find that the observed 
state is consistent with tidal theory and stellar evolution. 
For a star born with T e g < 6500 K, it could develop 
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a convective zone and lose angular momentum through 
"magnetic braking" . So even when it rotates quickly at 
birth, it will spin down quickly and thus have a longer pe- 
riod at ZAMS. While for a star born with T cS > 6500 K, 
it will be fully radiative and could not lose angular mo- 
mentum through "magnetic braking" . Thus it will rotate 
much more quickly. From the evolution track in FigJS]we 
could see that the effective temperature of TYC 2087- 
00255-1 is initially around the transition point ~ 6500 K, 
and its evolution is difficult to know a priori. If fully ra- 
diative, we expect it to rotate quickly while on the main 
sequence, while if convective, it would have spun down 
and will rotate more slowly. Here we consider both sce- 
narios and find that we cannot rule either out, at least 
using a simple, coupled model for the expansion of the 
primary and the tidal evolution. 

As star evolves, its radius expands and hence we ex- 
pect the star to spin down via conservation of angular 
momentum. If the rotational frequency of the star is cj, 
then its time rate of change due to expansion is 



(1) 



duj 2lu dR 

-jj-\exp = 

The rate at which ui changes is therefore encapsulated in 
dR/dt, which we can derive from stellar evolution mod- 
els. Using the "Y 2 " models (jDemarque et al.ll2004| ). we 
fit a third order polynomial to the radius as a function 
of time: 

= 1.0275 + 0.1661*' - 0.0619i' 2 + 0.01147i' 3 (2) 

Rq 

where t' is the age of the star in Gyr. Differentiating 
with respect to if we find 



dR* 
~dP~ 



= 0.1661 - 0.1238*' + 0.03441i /: 



(3) 



from which we can solve Eq. (fTJ) . 

The tidal evolution is considerably more compli- 
cated as it depends on many more parameters, as 
well as the tidal model employed. Here we use 
the "constant-phase-lag" (CPL) model as described in 
iFerraz-Mello et al.l (|2008t) . This widely- used model as- 
sumes a constant phase offset between the location of 
the companion and the tidal bulge. The magnitude of 
the phase lag is jq, where Q is the "tidal quality fac- 
tor." Different values of Q have been proposed to ex- 
plain tidal evol ution of different sys t ems ( 10 7 in iZahnl 
(fl989l); 10 5 in IMeibom fc Mathieul j200l) ;10 6 to 10 7 
in iSchlaufman et al.l (|2010l) ; 10 s to 10 9 in iPenev et all 
(|201lD ). The speed of the evolution is also a func- 
tion of the Love number of degree two, which is a 
measure of the height of the tidal bulge. Rather than 
reproduce the set of six coupled differential equations 
that comprises the CPL model, the reader is referred to 
Ferraz-Mello et al.l (120081). Other tidal mode ls exist (e.g. 



Hutill98lULeconte et al.ll2010t lHansenllMoT) that make 



qualitatively different assumptions. As we are only inter- 
ested in demonstrating that the observed configuration 
is consistent with tidal theory, we limit our scope to the 
CPL mod el. We use the nume rical methods outlined in 
App. E of lBarnes etaTl (|2012l ). 

We can use the best fit parameters in the CPL model, 
but we also must set Q and k 2 , and the moment of iner- 



tia constant, or "radius of gyration" r g . A wide range of 
values have been proposed for Q with 10 6 being a stan- 
dard choice (e.g. Uackson et aLll2008l l2009h . We set k 2 
to 0.5, which is arbitrary since the tidal evolution actu- 
ally depends on the quotient of Q and k 2 , and currently 
we cannot disentangle the two. We set R g to 0.35, con- 
sistent with theoretical e xpectations for solar-like stars 
(IClaret fc GimenezlH990l) . We assume the companion is 
tidally locked, use its minimum mass, and has a radius 
of 1 Rj U p. We set both bodies' obliquities to 0. With 
these choices, we may integrate the tidal evolution of the 
system forward, tracking the spins, obliquities, orbital 
period, and orbital eccentricity. 

Because the stellar rotational frequency depends on 
both its internal and tidal evolution, we include both 
effects in our model. Furthermore, we assume that both 
effects act independently, i.e. no feedbacks are present. 
This so rt of coupling has been applied to stellar binaries 
befor e (IZahn fc Boucheti 119891: iKhaliullin fc Khaliullinal 
120111: iGomez Maaueo Chew et al.ll2012[ ). but this may 
be the first time it has been done for a star with a BD 
companion. 

We first consider the convective case, i.e. slow ini- 
tial rotation. In Fig. [12j we show one plausible history 
for this system. We start the integration at an age of 
221 Myr when the star is on the main sequenc e, at which 
point the "Y 2 " model (|Demarque et al.ll2004D predicts a 
radius of 1.08 Rq. The orbit begins with a period of 
9.5 days, and an eccentricity of 0.24. The initial stel- 
lar rotation period is 40 days. We set the stellar Q to 
4 x 10 6 . Based on the radial evolution and the observed 
uncertainties in radius, we estimate the age of the sys- 
tem to be between 4 and 5.1 Gyr, and with a nominal 
age of 4.75 Gyr, corresponding to the best fit radius of 
1.64 Rq. This interval is shaded gray. This age estimate 
is different from that in §3.3[ which is estimated from 
T e s and log(g) on the HR diagram. We also considered 
a wider range of configurations and find that, for plau- 
sible primary spin periods, the primary's Q must lie in 
the range 3 x 10 6 - 6 x 10 6 in order to produce a system 
consistent with the observations. 

In the top panel, the evolution of the orbital period 
P or b is shown by the solid curve. The dashed curve is 
the best fit orbit, whose uncertainty is less than the curve 
thickness. The next panel down shows the eccentricity 
evolution. The line styles are the same as before, but 
now the la observational uncertainty in e is denoted by 
the horizontal dotted lines. Next is the stellar rotation 
period, shown in the same format as the previous panels. 
Fourth is the evolution of the stellar radius, as given by 
the fit in Eq. (|2|). Note that all these parameters pass 
through their best fit values at our nominal age estimate. 

The final panel shows the ratio of the stellar spin evo- 
lution from tides to that from radial expansion (Eq. pQ ) , 



dtu I 



W EE 



duj I 

dt \ ex p 



(4) 



where the denominator is the time rate of change of the 
rotational frequency due to tides. The bottom panel of 
Fig. [T2l shows the evolution of W. Although W evolves, 
it is always (1) negative, and (2) more negative than —1. 
These features indicate that the tidal torques oppose the 
radial expansion, and dominate. Therefore, the star is 
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Fig. 12. — Evolution of various properties of the MARVELS- 
4 system assuming a slow initial rotation period for the primary. 
Top: Orbital period. Top Middle: Eccentricity. Middle: Primary's 
rotation period. Bottom Middle: Primary Radius. Bottom: Ratio 
of the time rate of change of the primary's spin period from tides 
to that from expansion, c.f. Eq. Bt . 
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Fig. 13. — Evolution of various properties of the MARVELS-4 
system assuming a fast ini tial rotation period for the primary. The 
format is the same as Fig. 1121 

spinning up, and will continue to do so until it becomes 
tidally locked. (The discontinuity at 2.6 Gyr is due to 
passage through the 2 : 1 spin-orbit resonance.) 

We now turn to the possibility that the primary's rota- 
tion was initially fast. At first glance, it may appear that 



the observed system is inconsistent with such a history 
because the primary's rotation period is currently longer 
than the orbital period, implying the system has passed 
through the 1:1 spin-orbit resonance. In that case, the 
primary may have become tidally locked. However, if 
the expansion is rapid enough, and/or the tidal Q large 
enough, the primary could pass through this state and 
avoid permanent capture. In Fig. [T3l we show such a 
configuration. In this case, the primary's rotation period 
is initially 6 days, the eccentricity is 0.228, the primary's 
rotation period is 6 days, and its tidal Q is 3 x 10 7 . As 
before the model predicts a system that could evolve to 
the observed state. 

The evolution of W shows the complex evolution that 
may have occurred. Initially W > meaning that both 
tides and expansion act in the same direction, increasing 
the rotational period in this case. Very quickly the sys- 
tem passes through the 3:2 spin orbit resonance and the 
phase lags change, producing the sudden drop in W. For 
the next 2 Gyr W increases slightly because the orbit is 
shrinking and increasing the tidal torques. As the pri- 
mary approaches the 1:1 resonance, the torque decreases 
and W decreases accordingly. At 3.2 Gyr, the primary 
passes through the 1:1 resonance and now the tides act to 
speed up the rotation and hence W drops to less than 0. 
Right after the resonance crossing, the tidal torques are 
weaker, because the rotation is close to the equilibrium 
value (the orbital period). However, as the primary con- 
tinues to expand and slow down, the disparity increases 
and the tidal torques grow larger. Thus, \W\ increases 
because the tidal torque increases faster than radial ex- 
pansion slows the spin period. As W never becomes more 
negative than -1, the radial expansion dominates the evo- 
lution and rotational slowdown continues. 

We therefore have two competing, plausible evolution- 
ary models for this system. We have neglected some 
effects, such as magnetic breaking and coupling between 
the expansion and tidal evolution, and hence we do not 
express a preference for either model. However, given 
the uniqueness of this system (evolving F star, BD/low 
mass stellar companion, W ~ 0), it may be fertile ground 
for further exploration and insight into tidal processes on 
stars with little or no convective envelope. 

In § 14.41 we have found a point source near TYC 2087- 
00255-1 using AO imaging. Interactions between 
MARVELS-4b and the tertiary may serve the purpose 
of bringing MARVELS-4b from its birth place to a 
tight orbit in the early history of this system through 
Kozai-Lidov mechanism if the tertiary is indeed as- 
sociated with TYC 2087-00255-1 and if the initial 
mutual orbit inclination angle between M ARVELS-4b 
and the tert iary is 39.2° < S 23 < 141.8° ([Kozail Hl6l 
iLidovl I1962I ). This mechanism combined with tidal 
friction have been proposed to explain formation of close 
binaries in triple-star system ([Mazeh fc Shahaml 1 1979 ; 
Kiseleva et alJll998tlEggletoii fc Kiseleva-Eggletorj200li 
Egglc ton fc Kisseleva-Eggletonl 120061 : iTokovinin et al.l 
20061 : iFabrvckv fc Tremaind I2007D and formation of 
close-in Jupit er mas s exoplanets ("Hot Ju pi ters" ; 
Wu fc Murray! 120031 IFabrvckv fc Tremaind 120071 
Wuet al.1 120071 INaoz et all 1201 ID . After MARVELS-4b 
has been brought to a tight orbit, tidal force from 
the primary takes over and MARVELS-4b follows 
the evolution illustrated above in Fig. [12] and Fig. [13] 
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qualitatively. 

5.3. Expected Stellar RV Jitter 

Starspots and motions of the stellar surface are as- 
trophysical sources of noise that can interfere with 
searches for companion RV signals. These sources 
are commonly refe r red t o as "jitter" , first notice d 



1987) 



by IGunn fc Griffin] (|1979T > a nd ILupton et al 
and subsequently explored by ISaar fc Donahue! (11997). 
ISaar et all (IT998T). I Wrighfl (120051 ). ILagrange et all (|2009T) 
and llsaacson fc Fischer! (I2010D. Us i ng th e analyti- 
cal relation given in ISaar fc Donahue! (Jl997), we esti- 
mate the expected RV jitter for TYC 2087-00255-1 as 
6.5f 9 v ro t sin i, where / is the flux change in percent and 
v rot sin i is the projected rotational velocity. We have de- 
rived v ro t sin i from our spectra using two different meth- 
ods (see £|3.4[) . The final combined value for v rot smi is 
9.9 ± 0.8 km s _1 . The percentage change of the stellar 
flux is estimated to be / ~ 1.8 from the Super WASP 
photometry data, thus the expected RV jitter is ~ 106 
m s _1 . 

During the joint Keplerian orbital fit to the MAR- 
VEL+SARG RV data, a RV jitter - 112 m s" 1 was 
needed to account for the extra noise in our RV mea- 
surements, which matches the expected jitter. However, 
when we use the MARVELS RV data only, we yield a 
'jitter' term Ujitter = 152 m s _1 in our Keplerian orbital 
fit, which is bigger than the expected RV 'jitter' arising 
from the stellar activity. This implies that uncharacter- 
ized systematics remain in the MARVELS RV data, and 
may likely dominate the stellar RV 'jitter'. 

6. SUMMARY 

In a search through the first two years of SDSS-III 
MARVELS data, we discovered MARVELS-4b, a can- 
didate BD companion to the V ~ 10.6 star TYC 2087- 
00255-1 with a velocity semi-amplitude of K — 3.571 ± 
0.041km s _1 and a short orbital period of 9.0090 ± 
0.0004d, yet with an eccentricity e = 0.226 ±0.011. Ad- 
ditional RV data from SARG observations confirm the 
doppler variability. High-resolution spectroscopic obser- 
vations indicate that the host star is a slightly evolved 
subgiant with T off = 5903 ± 42 K, log^ = 4.07 ± 0.16, 
and [Fe/H] =-0.23 ±0.04, with an inferred mass of M* = 
1.16 ± 0.11 M Q . The minimum mass of MARVELS-4b is 
40.0 ±2.5 Mj up , implying that it is most likely in the BD 
regime. A 13.16 day periodic signal is found in the Super- 
WASP photometry data, which is likely due to rotational 
modulation of starspots on the host star, and indicates 
that this star-BD system is not tidally synchronized. A 
second possible companion is found 643 ±10 mas away 
from TYC 2087-00255-1 using K' AO imaging. Its asso- 
ciation with the primary star could be verified by future 
proper motion measurements. Ca II H and K core emis- 
sion indicates that the subgiant is chromospherically ac- 
tive at a level unusual for subgiants. Tidal interactions 
between the star and BD could have spun up the star 
and make it active. 
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